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ABSTRACT 

Current models of galaxy formation lack an efficient and physically constrained mech- 
anism to regulate star formation (SF) in low and intermediate mass galaxies. We argue 
that the missing ingredient could be the effect of photoionization by local sources on 
the gas cooling. We show that the soft X-ray and EUV flux generated by SF is able to 
efficiently remove the main coolants (e.g., He + , 4+ and Fe 8+ ) from the halo gas via 
direct photoionization. As a consequence, the cooling and accretion time of the gas 
surrounding star-forming galaxies may increase by one or two orders of magnitude. 
For a given halo mass and redshift, the effect is directly related to the value of the star 
formation rate (SFR). Our results suggest the existence of a critical SFR above which 
"cold" mode accretion is stopped, even for haloes with M V i r well below the critical 
shock-heating mass suggested by previous studies. The evolution of the critical SFR 
with redshift, for a given halo mass, resembles the respective steep evolution of the 
observed SFR for z < 1. This suggests that photoionization by local sources would 
be able to regulate gas accretion and star formation, without the need for additional, 
strong feedback processes. 
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1 INTRODUCTION 

Despite many successes, current models of galaxy formation 
have still several difficulties in reproducing the galaxies as 
we see them today. The gas "overcooling" in cosmological 
simulations is considered as one of the main sources of this 
discrepancy with observations. Feedback mechanisms like 
supernovae (SN) or galactic winds heating have been ad- 
vocated in order to suppress gas cooling and star formation 
(Dekel & Silk 1986). However, numerical models including 
these feedback effects have not been able, so far, to obtain 
realistic galaxies and to reproduce the steep decline of the 
star formation rate at z < 1 (see e.g., Crain et al. 2009). 

While this may be related to numerical limitations (that 
force the introduction of ad hoc prescriptions), it is worth- 
while to re-consider on a physical basis the gas cooling mech- 
anism and how the galaxies get their gas. As originally 
pointed out by Silk (1977) and Rees & Ostriker (1977), 
galaxy-formation is regulated by the relative value of the 
cooling and collapse timescales. Gas in haloes with masses 
below 10 M©, cools so efficiently that it always accretes, 
roughly, at the free-fall rate. This "cold" mode accretion 
(as opposite to the "hot" mode, where gas is stably shock- 
heated to the halo virial temperature) has been confirmed re- 
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cently by more sophisticated analytical and numerical mod- 
els (Dekel & Birnboim 2006, Ocvirk et al. 2008 and ref- 
erences therein). The physics of gas cooling is determined 
by the abundance of the ions that act as coolants for the 
free electrons in the gas. Thus, every process that changes 
these abundances will act directly on the galaxy-formation 
timescales. This simple idea has been the base of the seminal 
work by Efstathiou (1992), where it is clearly illustrated how 
the cosmic photoionizing background suppresses the cooling 
due to H° and He + for low density gas (see also Wiersma et 
al. 2009) . This effect is important to suppress the formation 
of (dwarf) galaxies within dark-matter haloes less massive 
than 1O 9 M0 but it does not substantially affect cooling and 
accretion of virialized, dense gas in more massive haloes (in- 
deed, the majority of numerical models already include this 
effect in the cooling function). 

In this Letter, we consider an ingredient neglected so 
far: the effect of the local photoionizing flux from star for- 
mation on gas cooling and accretion. Star-forming galaxies 
are conspicuous sources of soft X-ray and EUV photons. In 
Section 2, we discuss their spectral energy distribution. In 
Sections 3 and 4, we calculate the effect of the local pho- 
toionization flux on gas cooling and galaxy formation. We 
conclude in Section 5. We use a flat ACDM cosmology with 
h = 0.7, On = 0.27, Q.A = 0.73, Q b = 0.046, X He = 0.25, 
and we denote with "pKpc" distances in physical kpc. 
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Figure 1. Net cooling rates (defined as the absolute value of the radiative cooling minus the photo-heating rate) as a function of gas 
temperature for three different metallicities (see labels in panels). The gas density is rin = 1.6 X 10~' 5 cm " 3 (A ~ 10 3 at z = 1). The 
light solid lines show the cooling rate in the Collisional Ionization Equilibrium (CIE) case, when no photoionizing flux is present. The 
dark solid lines represent the net cooling rate (note that heating dominates over cooling at low temperatures) in the presence of the 
cosmic UV background (HM05) alone at z = 1, while the other lines show the effect of adding the local photoionizing flux from star 
formation with log(C7*) = [—1,0, 1,2] (see labels). These values can be directly translated to SFRs using Eq[T]or Eq[2] 



2 THE PHOTOIONIZING FLUX OF 
STAR-FORMING GALAXIES 

The onset of star formation episodes within galaxies pro- 
duces a large amount of ionizing photons associated with a 
range of different processes. The EUV part of the spectrum 
is dominated by the continuum of massive, hot and young 
stars. Conspicuous X-ray emission is produced by the me- 
chanical energy released into the Interstellar medium (ISM) 
during the starburst (Heckman et al. 1995, Strickland et al. 
2004), by non-thermal processes associated with supernovae 
explosions, and by massive X-ray binaries (Grimm et al. 
2003, Persic et al. 2004). In particular, the soft X-ray flux 
correlates well enough with the star formation rate (SFR) 
obtained from other indicators (e.g., the far infrared lumi- 
nosity) to be used as a proxy to derive the SFR for actively 
star- forming galaxies (e.g., Rosa-Gonzalez et al. 2009 and 
references therein). 

For this reason, we use the spectral energy distribution 
(SED) derived from models that are calibrated empirically 
to reproduce the relation between SFR and soft X-ray emis- 
sion. In particular, we use the model of Cervino et al. 2002 
(http://www.laeff.inta.es/users/mcs/SED), that is consis- 
tent with observations assuming: i) relatively young star for- 
mation episodes (age of about 5 Myrs) and, ii) a few per cent 
efficiency (~ 5%) in converting mechanical energy from the 
starburst into ISM heating and X-ray emission (Mas-Hesse 
et al. 2008). Assuming a Salpeter IMF and Z = Z , the 
relation between the SFR and the soft X-ray (0.5-2.0KeV) 
emission for this model is: L soft x = 3 x 10 40 [SFR/(1 M 
yr -1 )] erg s _1 . We use this relation to calibrate the overall 
normalization of the intrinsic SED as a function of the SFR. 

For typical star-forming galaxies, the soft X-rays flux is 
not substantially affected by ISM absorption. Thus, our re- 
sults would be independent on absorption effects if the main 
coolants are ions with ionization potential corresponding to 
the energy of soft X-rays photons. As we show in Section 
3, this is indeed the case for gas with T ^ 10 5 K. For lower 
temperatures, H° absorption becomes more important and 
we must consider how many photons with frequency around 
the Lyman Limit can actually escape the ISM. In particular, 



we use a frequency-dependent escape fraction {f" sc ) of the 
form/e SC = [/^ + (l-/«c)«!p(-7V')]) where f££ is the ab- 
solute escape fraction at the Lyman Limit, t„ = aviV(H ) is 
the neutral hydrogen optical depth and ov the correspond- 
ing cross section. This is compatible with a leaky-absorbers 
model (see, e.g. Ferland & Mushotzky 1982) for the ISM 
with a covering factor of (1 — fesc)- We fix = 0.05, 
compatibly with observational estimates for local starburst 
galaxies (e.g., Bergvall et al. 2006, but see Grimes et al. 
2009) and high-redshift LBGs (Steidel et al. 2001, Shapley 
et al. 2006). The value of A^H ) determines the hardening of 
the spectrum around the Lyman Limit. We anticipate that 
this parameter has a little effect on our results and we fix 
N(B°) = 10 20 cm~ 2 . Note that He and dust may be also 
an important sources of opacity. We found that these effects 
become significant in reducing the EUV flux for gas column 
densities higher than iV(H) = 10 21 cm -2 , assuming Galactic 
ISM abundances. 

The relation between the SFR and the total luminosity 
of ionizing photons (y > ^ll) for the assumed SED is given 
by: L ion = 3 x 10 41 [SFR/lMoyr -1 ] erg s" 1 . We stress that, 
even if Li 0n depends on /i^, the overall effect of the ionizing 
photons on gas with T > 10 s K is almost independent on 
this parameteiQ. Nevertheless, we use Li OU (in terms of the 
ionization parameteiQ) as a convenient way to parametrise 
our results in the rest of this Letter. 



3 LOCAL PHOTOIONIZING FLUX AND GAS 
COOLING 

We calculate now the effect of the ionizing photons on the 
cooling of a gas patch in the galaxy halo. For the assumed 
SED (Section 2), the ionization parameter is given by: 

For instance, we have verified that a much lower escape fraction, 
fesi = 0.005, would change the cooling rates at T ^ 10 5 K only 
by few percents with respect to the values presented in Section 3. 
2 defined as the ratio between the density of photons with v > 
I'LL and the hydrogen density at a given distance from the source. 
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U, = 0.03 x 



SFR 



lM yr- 



7 pKpc 



10 3 



1+2 



(1) 



where r and A are, respectively, the distance from the galaxy 
and the gas overdensity. For comparison, the ionization pa- 
rameter given by the cosmic UV background (Haardt & 
Madau 2005, in preparation; HM05) at z = 1 is Uhm = 
7.5 x 10~ 3 for A = 10 3 . Note, however, that the pho- 
ton spectral distribution of the HM05 to the SED of the 
star- forming galaxies are substantially different, especially 
around the peak of quasar activity (z ~ 2). Using the halo 
virial relations (see e.g., Mo & White 2002 and DB06 for 
an useful compilation) and an isothermal gas profile (i.e., 
A(r) tx r -2 ), Eq[T] translates into: 

SFR 



[/* = 0.03 x 



lM yr- 



(MiT) 



2/3 



0.05 



h 



l + z 



(2) 



where M\{ T is the virial mass in units of 10 11 



M and fb is 

the gas fraction inside the halo. Note that, for an isothermal 
profile, Ut is independent on the distance from the galaxy. 

How does the inclusion of the local photoionizing flux 
change the gas cooling? For a gas with primordial composi- 
tion, we could estimate the effect on the two main coolants, 
H° and He + (mainly via HI and Hell Lya emission), di- 
rectly from the value of the ionization parameter. However, 
in a more realistic situation, traces of metals are present 
and dominate the cooling at T 5 > 1 (where T 5 = T/[W 5 K]) 
via collisional excitations of several ions. In order to get the 
detailed ionic balance in the presence of an EUV/X-ray ra- 
diation field and, ultimately, the heating and cooling rates, 
we use the publicly available package Cloudy (version 08.00 
of the code last described by Ferland et al. 1998), which 
contains most of the atomic processes that are important 
for our temperature range. 

The results of the photoionization model are presented 
in Fig.l, where we show the absolute value of net cool- 
ing/heating rate, A nei = A — H (A and H are, respectively, 
the radiative cooling and the photo-heating rate in erg cm -3 
s" 1 ) as a function of the gas temperature for three differ- 
ent metallicities. In this plot, the gas density is fixed at 
tlh = 1-6 x 10 -3 cm -3 , corresponding to A = 10 3 at z = 1. 
In each panel, the thick solid line represents the net cool- 
ing rate (note that heating dominates over cooling at low 
temperatures) in the presence of the HM05 alone, while the 
other thick lines show the effect of adding the local photoion- 
izing flux from star formation with log(£/*) = [—1,0,1,2] 
(see labels). These values correspond roughly to SFR ~ 
[3, 30, 300, 3000] for r ~ 7 pKpc (EqQJ or Mjf" ~ 1 (Eq0. 

In the presence of a photoionization flux, the main cool- 
ing peaks are due to He + (at T5 ~ 1), 4+ (at T5 ~ 2, via 
O V 630A line), Ne 5+ (at T 5 ~ 5, via Ne VI 400A line) and 
Fe s+ (at T 5 ~ 10, via Fe IX 169A line). Clearly, the local 
photoionizing flux strongly reduces the net cooling rate at 
T 5 ~ 1 for U, ^ 1 and at T 5 ~ 10 for U* ^ 10. This is 
due both to the suppression of coolant ions via photoioniza- 
tion and the injection of energy following the photoioniza- 
tion process (i.e., photoheating). In the extreme case of a 
highly overionized medium (i.e., U* = 100, dot-dashed lines 
in Fig.l) all the coolant ions are removed and the main cool- 
ing process if free-free emission, which is highly inefficient at 
low temperatures, but is important for T > lO 7 ^ gas (e.g., 
the Intra Cluster Medium). 
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Figure 2. Cooling time as a function of gas temperature (or, 
equivalently, M v i r ) in the presence of the cosmic UV background 
only (SFR= 0) and for different values of the SFR (in units of 
Mq yr -1 ) and redshifts (see labels). The gas density is calculated 
(at 0.1-R„i r ) using the halo virial relations taking into account the 
redshift evolution and the mass dependence of the concentration 
parameter (see text for details) . The gas metallicity is fixed at Z = 
O.O3Z0. The horizontal lines give the value of the compression 



time, t co „ 



defining the transition between "cold" and "hot" 



accretion mode at the corresponding redshift. 



4 IMPLICATION FOR GALAXY FORMATION 

The cooling time of a patch of gas with density nn and 
temperature T is defined by: 



tc 



2 A net (T,n H , Z, U*,z) 



(3) 
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We compare t coo i to the compression time of the gas as it 
falls within the dark-matter halo (tcomp, roughly equal to 
few times the gas free fall time ), following Dekel & Birn- 
boim 2006 (DB06). Note that the last equation is valid in 
the innermost part of the halo (at about 10% of the virial 
radius). As extensively discussed by DB06, the relation be- 
tween these two time-scales around r = O.litw seems to be 
the origin of the transition between "cold" and "hot" mode 
accretion. The actual value of the compression time depends 
on the details of the shock (e.g., the shock speed) and the gas 
profile. Nevertheless, the simple condition t coo i > t cornp (as 
defined above) seems to produce very similar results to the 
one obtained from numerical simulations (e.g., Ocvirk et al. 
2008), where no assumptions about the shock properties are 
made, providing that a low gas metallicity (Z ~ 0.03Zq) 
is assumed. Therefore, in the rest of this Letter, we use 
Z = O.O3Z and r = 0.1R vir . 

In Fig. 2, we show the gas cooling time as a function 
of the virial shock temperature (or virial mass) for three 
different redshifts (see labels in panels). In each panel, the 
black solid line represents the cooling time in the presence 
of the HM05 only, while the horizontal lines show the value 
of t comp at the given redshift. We have computed the gas 
overdensity at r = 0.1-R^j,- using the halo virial relations, 
assuming an isothermal profile for the gas, a halo gas frac- 
tion fb — 0.05 (consistently with DB06), and taking into 
account the evolution of the haloes concentration with red- 
shift and mass (using the fit of Bullock et al. 2001). Similarly 
to DB06, where local sources are not included, the transition 
between "hot" and "cold" mode happens at the critical halo 
mass Merit ~3x lO n M , for Z = 0.03Z Q , independently 
of redshift. The other lines in Fig. 2 show the cooling time 
in the presence of the local ionizing flux for different values 
of the SFR in units of M@ yr _1 (see labels). It is clear that 
(especially at low redshifts) the modality of gas accretion de- 
pends on the presence of star-formation activity: at a given 
redshift, M cr it is actually determined by the value of the 
SFR. Note that, at the lower mass end (Mfr < 0.5), cool- 
ing is completely suppressed (and the gas is actually slightly 
heated above T v i r ). 

In Fig. 3, we show the redshift evolution of the critical 
SFR (i.e., the SFR at which t coo i — t comp for r = 0.1R vir and 
"cold" mode accretion is stopped; solid red lines) for differ- 
ent halo masses. Note that the critical SFR increases with 
halo masses for Mil r < 1 (upper panel) while the trends 
is inverted for M^{ r > 1 (lower panel). This is simply due 
to the shape of the cooling function or, equivalently, to the 
relation between the cooling time and halo mass (see Fig- 
ure 2). For illustrative purposes, we overlay the observed 
SFRs of galaxies with total stellar mass of the order of 
M, ~ 2 x 1O 9 M (upper panel) and M, ~ 1O 1O M (lower 
panel) from various sources (see caption). To facilitate the 
comparison, we converted the stellar masses to an approxi- 
mate M v ir range extrapolating the results of Mandelbaum et 
al. (2006). Note that large uncertainties are associated with 
this conversion, especially for the low mass bin. Nevertheless, 
it appears remarkable that the critical SFR (that basically 
depends only on the cooling physics) is close to the observed 
value. Also, the slope of the redshift evolution appears very 
similar to the observed relation. Note that the critical SFR 
represents the instantaneous value of the SFR that is needed 
in order to stop "cold" mode accretion, while the observed 
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Figure 3. Evolution of the critical SFR (the SFR at which 
tcooi = tcomp for r = 0.1R v ir, and "cold" mode is stopped) for 
different halo masses (solid red contour lines). The gas density is 
calculated using the halo virial relations taking into account the 
redshift evolution and the mass dependence of the concentration 
parameter (see text for details). The gas metallicity is fixed at 
Z = O.O3Z0. In the upper (lower) panel, contours are equally sep- 
arated by 0.1 (0.2) in M^ r . Note that the critical SFR increases 
with halo masses for M^"~ < 1 (upper panel) while the trends is 
inverted for M?^ r > 1 (lower panel), as expected given the rela- 
tion between cooling time and halo masses shown in Fig. 2. For 
illustrative purposes, we overlay the observed SFRs of galaxies 
with M, ~ 2 X 10 9 M Q (upper panel) and A/* ~ 10 10 M Q (lower 
panel) obtained by the SDSS survey (crosses), Zengh et al. 2007 
(open circles) and Damen et al. 2009 (circles). The observed stel- 
lar masses have been converted to an approximate M v i r range 
extrapolating the results of Mandelbaum et al.(2006). 

values represent the average SFRs of the galaxies (that may 
also be at different stages of evolution). In this respect, the 
critical SFR should be regarded as a sort of maximum and 
for this reason we might expect the observed values to be 
below it, as in Fig. 3. 



5 DISCUSSION AND CONCLUSIONS 

We have argued that the local photoionizing flux produced 
by star-forming galaxies could play a fundamental role for 
the cooling and accretion of the halo gas. In particular, we 
have found that the soft X-rays and EUV photons pro- 
duced by star formation activity, are able to remove from 
the gas important coolants (e.g., 4+ and Fe 8+ ) via direct 
photoionization. The reduced cooling rate (and, partially, 
the increased photoheating) result in a much longer cooling 
time for the accreting gas. As a consequence, "cold" mode 
accretion may be stopped for haloes well below the critical 
halo mass for shock-heating, M cr it ~ 3 x lO n M0, inferred 
neglecting local photoionization sources. 

Specifically, we have found that (for a given halo mass 
and redshift) there is a critical SFR corresponding to the 
transition between "cold" and "hot" mode accretion. The 
evolution of the critical SFR with redshift, for a given halo 
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mass, is of the form (1 + z) a with a ~ 3, resembling 
the corresponding steep evolution of the observed SFR for 
z < 1 (e.g., Damen et al. 2009). This intriguing result sug- 
gests that, at least for this redshift range and for haloes 
with 10 9 - 5 < (M vir /M ) < 10 11,5 , photoionization by local 
sources may represent an important mechanism to regulate 
gas accretion and star formation helping to relieve the need 
of additional, strong feedback processes (e.g., galactic winds 
or mechanical SN feedback) . This may help to solve a long- 
standing problem of numerical models of galaxy formation: 
the effective quenching of star formation in intermediate and 
low-mass haloes. 

How sensitive are these results to variation in metallic- 
ity and densities assumed here? Following DB06, we used a 
simple model based on the halo virial relations and the evo- 
lution of the concentration parameter to derive the densities 
that enter into the definition of the cooling time. This ap- 
proach gives very similar transition masses between "cold" 
and "hot" mode accretion with respect to hydro-dynamical 
simulations (e.g., Ocvirk et al. 2008). These models also find 
that "cold" mode accretion is typically associated with low 
metallicity gas, as we have assumed. Nevertheless, single ob- 
jects may have different densities and metallicities with re- 
spect to the "average" values used here. In order to asses 
what is the effect of an enhanced or decreased density on 
our results, we have rescaled the densities given by the halo 
virial relations by a factor f e . We found that the critical 
SFR scales roughly like / e 2 (for MJf - 1 and 0.5 < f e < 5), 
similarly to what expected in the CIE case. In the case of a 
more metal enriched gas the critical SFR also increases, up 
to a factor of two for Z — O.IZq. These effects may introduce 
a scatter into the relation between the critical SFR and halo 
masses (changing the normalization of the curves presented 
in Figure 3), however they do not change its steep evolution 
with redshift, that is in good agreement with observations. 

What is the effect of additional sources of photoion- 
ization and heating? AGN activity is sometimes connected 
with star-forming galaxies, boosting by few order of mag- 
nitudes the soft-X ray emission (see e.g., Alexander et al. 
2005, Georgakakis et al. 2007). In this respect, the results 
presented in this Letter should be regarded as the minimum 
effect due to star-forming galaxies when the AGN are not 
present. Moreover, AGN may be also able to strongly re- 
duce cooling in the "hot" mode. For a tipical luminosity of 
isoftx = 10 43 erg s _1 the corresponding ionization parame- 
ter is high enough to significantly affect the cooling rate also 
for T > 10 6 K (see lines with £7* > 1 in Fig. 1). This is con- 
sistent, e.g., with the findings of Oh (2004), that discussed 
this effect as a potential mechanism to suppress Fe emission 
lines from cooling-flow clusters. 

Star formation should be followed by the release of large 
amount of energy via SN explosions or stellar winds (Dekel & 
Silk 1986) that we have not included (as a heating source) 
in our calculations. However, the effective amount of en- 
ergy released into the halo gas and the physical details of 
the process are affected by large uncertainties connected, 
e.g., with the mass loading factor of the galactic winds. We 
stress that the effect discussed in this Letter is substan- 
tially different from SN and winds feedback, even if they 
are all connected with star-formation activity: SN feedback 
heats the gas while photoionization from local sources re- 
duces the cooling rate removing coolant agents (and also 



adds energy through photo-heating). These two effects may 
work together, reducing the critical SFR or e.g., helping the 
re-heated gas from SN feedback to remain hot preventing it 
to cool and to fall back onto the galaxy. 

Despite the limitations of the simple analytical argu- 
ments used here, photoionization by local sources certainly 
represents a fundamental ingredient that should be added 
into current galaxy-formation models. In a companion pa- 
per, we will show how this effect changes the galaxy proper- 
ties with the help of cosmological, hydro-dynamical simula- 
tions. For the time being, the simple estimates presented in 
this Letter may represent a further step towards a better un- 
derstanding of the processes that regulate galaxy-formation 
and produces the galaxies as we see them today. 
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